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ABSTRACT 

The location of the Scutum Red-Supergiant (RSG) clusters at the end of the Galactic Bar makes 
them an excellent probe of the Galaxy's secular evolution; while the clusters themselves are ideal 
testbeds in which to study the predictions of stellar evolutionary theory. To this end, we present a 
study of the RSGs' surface abundances using a combination of high-resolution Keck/NIRSPEC H-band 
spectroscopy and spectral synthesis analysis. We provide abundance measurements for elements C, O, 
Si, Mg, Ti, and Fe. We find that the surface abundances of the stars studied are consistent with CNO 
burning and deep, rotationally enhanced mixing. The average a/Fe ratios of the clusters are solar, 
consistent with a thin-disk population. However, we find significantly sub-solar Fc/H ratios for each 
cluster, a result which strongly contradicts a simple extrapolation of the Galactic metallicity gradient 
to lower Galacto-centric distances. We suggest that a simple one-dimensional parameterization of 
the Galaxy's abundance patterns is insufficient at low Galactocentric distances, as large azimuthal 
variations may be present. Indeed, we show that the abundances of O, Si and Mg are consistent with 
independent measurements of objects in similar locations in the Galaxy. In combining our results 
with other data in the literature, we present evidence for large-scale (~kpc) azimuthal variations 
in abundances at Galacto-centric distances of 3-5 kpc. While we cannot rule-out that this observed 
behaviour is due to systematic offsets between different measurement techniques, we do find evidence 
for similar behaviour in a study of the barred-spiral galaxy NGC 4736 which uses homogeneous 
methodology. We suggest that these azimuthal abundance variations could result from the intense 
but patchy star formation driven by the potential of the central bar. 

Subject headings: open clusters & associations, supergiants, stars:evolution, stars:late-type, Galaxy: 
abundances, Galaxy: evolution, Galaxy: disk 



1. INTRODUCTION 

Massive young clusters are powerful natural laborato- 
ries with which to study many different aspects of as- 
trophysics. They contain large numbers of coeval stars 
with similar initial chemical compositions. Their youth 
and large initial masses mean that they contain many 
massive stars, and hence are intrinsically very bright ob- 
jects. Thus, such clusters can be used to make critical 
tests of stellar evolution, while they also trace recent star- 
forming history and local abundance patterns. 

Two such clusters exist in the constellation of Scutum, 
at distances of ^6 kpc. They are among the most mas- 
sive young clusters in the Galaxy, with initial masses 
in excess of 10 4 M Q . Their ages are tuned in such a 
way that their most massive remaining stars are in the 
Red Supergiant (RSG) phase, and together these clusters 
contain ^20% of all known RSGs in the Galaxy - clus- 
ter 1 (hereafter RSGC1) contains 14 RSGs , while cluster 
2 (RSGC2) contains 26 (jFiger et alj 120061: iDavies et all 
[2lm[2008aL hereafter F06, D07, D08). 

The clusters are located in the Galactic Plane, and 
have high visual extinctions (Ay « 25 and 16 for RSGC1 
and RSGC2 respectively). However, the fluxes of the 
RSGs peak at ~lyLtm, making them extremely bright 
in the near-IR where the extinction is reduced. Hence, 



the RSGs themselves make excellent tools with which to 
probe the global properties of their host clusters. 

The clusters' distances from Earth and Galactic longi- 
tude place them at Galacto-centric distances of ^3.5 kpc, 
and at an azimuthal angle 4>gc = 45 ± 10° relative to 
the Sun-Galactic Centre axis. This pla ces them at the 
tip of the Galactic Bar as described by iBeniamin et al.l 
(2005), which has half-length 4.5kpc and 4> GC = 44±10°. 
This raises the intriguing possibility that the RSG clus- 
ters were formed in a region- wide starburst event, trig- 
gered by the interaction between the Bar and the disk. 
Aside from the two clusters, there is certainly plenty 
of other evidence for recent star-forming activity along 
this line-of-sight: D07 found several other red stars in 
the field of RSGC2, which while not belonging to the 
cluster, were likely to be supergiants and hence massive 
stars. This is consistent with the large number of RSGs 
found at the base of the Scu tum spiral arm found by 
lLopez-Corredoira et alJ (|1999f ). Also noted in D07 were 
other objects indicative of recent star-forming activity, 
such as the candidate supernova remnant IRAS 18369- 
0557 and the candidate massive evolved s tar IRAS 18367- 
0556; while iGotthelf fe Habernl (|2008f) have identified 
the TeV gamma-ray source HESS J1837-069 with pos- 
sibly two young pulsar wind-nebulae. 



2 



Davies et al. 



In this paper we assess the chemical abundances of the 
two clusters, using high-resolution iJ-band spectroscopy 
of a sample of RSGs in each cluster, in combination with 
LTE model atmospheres. Our method allows us to derive 
abundances for C, Fe, and the a-elements O, Ca, Si, Mg, 
and Ti. We use these chemical abundances to address 
two distinct topics. 

Firstly, we use the relative surface abundances to test 
the predictions of stellar evolutionary models. In massive 
stars, the primary route of H fusion is via the CNO-cycle. 
Due to the relative reaction rates involved in this cycle, a 
depletion of carbon occurs at the expense of nitrogen. As 
stars enter the RSG phase, their deep convective layers 
are expected to bring the products of nuclear burning 
to the surface. Thus, the relative C abundances in the 
atmospheres of the RSGs allow us to make critical tests 
of the latest stellar ev olution codes, such as th ose which 
include rotation (e.g. iMevnet fe Maeden 20001) and the 
inter action with binary companions (e.g. Eldridgc et al. 
HH). 

Secondly, we will use the relative abundances of Fe 
and several a-elements in the RSG clusters, in combina- 
tion with the clusters' coincidence with the Galactic Bar, 
to investigate the star-formation history of the Galaxy. 
The dominant source of a enrichment in the Interstellar 
Medium (ISM) is through type-II supernovae (SNe), i.e. 
the core-collapse of massive stars. On the other hand, 
Fe-peak elements are produced in type-la SNe - ther- 
monuclear explosions of low mass stars which reach the 
Chandrasekhar limit through accretion from a compan- 
ion. Hence, a-elements are enriched on short (~Myr) 
timescales, while Fe-peak abundances are increased over 
much longer timescales (~Gyr). 

In our Galaxy, we see variations in the [a/Fe] ratio: 
in the central Galactic Bulge, at Galacto-ccntric dis- 
tances of i?cc <^2kpc, the [a/Fe] ratio is found to be 
super-solar from analyses of Red Giants (iRich fe Origlial 
20M ICunha fe Smithll20M iFulbright et al.ll2006l . (20071 
Lecureur et al. 1 120071) . In the 'thin' Galactic disk at 



Rgc ^ 5kpc and scale height \h\ <;100pc, the ratio 
q/Fe is found to be closer to solar (|Bensbv et al.ll200~l 
iLuck et ai1l2006t and references therein). In addition, 
the relative abundances in the thin disk are a function of 
i?GC : there is evidence to suggest heavy elemental abun- 
dance levels increase at lower Galacto-centric distances, 
while the ratio [ a/Fe] increases i n the outer disk (e.g 
ILuck et al.ll200l lYong et al.ll2006h . 

These abundance patterns are commonly explained as 
being due to the distinct star-formation histories of the 
different Galactic environments. The Bulge was formed 
rapidly several Gyrs ago, in a starburst event which was 
too brief for the chemical evolution to be affected by 
Type-la SNe. In the disk, star-formation has continued 
over the lifetime of the Galaxy, allowing Type-la SNe to 
contribute to the chemical enrichment of the ISM. 

In the inner disk (Rqc Ss 5kpc), the abundance pat- 
terns are thought to be strongly influenced by the Bar. 
External galaxies with bars are found to have signifi- 
cantly flatter abundance gradients than t hose without 
(jMartin fe Rovlll994t IZaritskv et aflll994D . This is ex- 
plained as being due to radial motions induced by th e bar 
potential (|Friedli et alJ 11994 iFriedli fe Bend Il995f ). In 
addition, pile-up of material at the bar's resonances can 



create star-forming hot-spots which result in localised 
chemical enrichment. The RSG clusters themselves are 
evidence of such a localised starburst in the inner regions 
of our own Galaxy, and have Galactocentric distances 
similar to co-rotation (~3.5kpc) and the Outer Lindblad 
Resonance (^5kpc). 

In this paper, we will explore the chemical abundances 
of the inner disk from analyses of the RSGs in the two 
Scutum clusters. By comparing to other abundance mea- 
surements in the inner disk, we will probe the star form- 
ing history in a region which is pivotal in the Galaxy's 
secular evolution. 

The paper is organized as follows: in Sect. [2] we de- 
scribe our observations, data-reduction steps and anal- 
ysis techniques; in Sect. [3] we present our results and 
discuss the abundance patterns of the two clusters. We 
make quantitative comparisons of these results with stel- 
lar evolutionary models in Sect. [4j while in Sect. [5] we 
discuss the abundances of the two clusters in the context 
of those of the inner Galaxy. We conclude in Sect. [6l 

2. OBSERVATIONS, DATA REDUCTION & ANALYSIS 

Observations were taken during the night of 13th Aug 
2006 using NIRSPEC, the cr oss-dispersed echelle spec- 
trograph mounted on Keck-II (|McLean et al.l l"l995). The 
instrument was used in high-resolution mode with the 
0.576" x 12" slit and the NIRSPEC-5 filter. We used a 
dispersion angle of 63.0° and a cross-dispersion angle of 
36.72°. This achieved a spectral resolution of ^17,000 at 
select regions within the wavelength range of 1.5-1.7/zm. 

We used integration times of 20s, and observed each 
star twice at different positions nodded along the slit. 
We observed all 14 RSGs in RSGC1, and 12 RSGs from 
RSGC2. The ide ntifications o f the speci fic stars ob- 
served , taken from iFiger et all (|2006l ) and iDavies et al.1 
(pool , can be found in Table[3] In addition to the cluster 
stars, we also observed HD 171305 (spectral type B0 V) 
as a telluric standard. Flat-field images were taken with 
a continuum lamp, while for wavelength calibration we 
observed Ar, Ne, Xe and Kr arc lamps, to provide a large 
number of template lines in the narrow wavelength range 
of each spectral order. 

Removal of sky emission, dark current and bias offset 
was done by subtracting nod-pairs, and images were flat- 
fielded with the continuum-lamp exposures. The spectral 
traces produced by NIRSPEC are warped in both the 
spatial and dispersion directions, so before extracting 
the spectra from the reduced science frames each two- 
dimensional order was rectified onto a linear grid. The 
spatial warping was characterized by the spectral traces 
of the two stellar traces in a nod pair. The distortion 
in the dispersion direction was fitted using the arc- and 
etalon-li nes. For technical details of the rectification pro- 
cess, see IFiger et all (|2003t) . 

As we know the wavelengths of the arc-lines, rectifica- 
tion also wavelength calibrates the data. By rectifying 
the arc frames and measuring the residuals of the arc- 
line wavelengths, we found the wavelength solution to be 
accurate to ±4 km s _1 across all orders. 

Spectra were extracted by summing each channel in 
the rectified orders within ±5cr of the centre of each spec- 
trum. Cosmic-ray hits and bad pixels were identified by 
comparing the two nod spectra of each star, and were 
replaced by the median value of the closest 4 pixels. 
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We removed the H and He I absorption features of the 
telluric standard via linear interpolation either side of the 
line. The atmospheric absorption features in the science 
frames were then removed by dividing through by the tel- 
luric standard. Finally, the spectra were normalised by 
dividing through by the median continuum value. From 
featureless continuum regions in the final spectra, we es- 
timate the signal-to-noise to be better than 100 for all 
spectra. 

2.1. Data analysis 

In order to obtain chemical abundances from the spec- 
tra, we have developed a technique of generating syn- 
thetic spectra from cool star model atmospheres, and 
adjusting the input relative abundances to match the 
equivalent widths of cer tain diagnosti c lines . The code 
used is that presented inlOriglia et al.l (|1993l ). which was 



updated in lOriglia et all (|2002t 12003 1. The code has 



been used successfully to obtain abundances of Bulge 
giants (e.g. iRich fe Orighalfe OOS). young clusters domi- 
nated by RSGs (e.g. lLarsen et alJ 12008). and most re- 
cently to study t he chemical abundan des of RSGs in the 
Galactic Centre (|Davies et al.ll2008bfl . 

The code and methodology is explained in detail in 
the papers listed above. Briefly, the code uses the LTE 
ap proximation a n d mo del atmospheres from the grid 
of Uohnson et alJ |l980) which incorporate the effect of 
molecular blanketing. Thousands of near-IR molecular 
ro-vibrational and atomic transitions are included. Ini- 
tial estimates of the stars' luminosities, temperatures and 
surface gravities were estimated based on the results of 
D07 and D08, while the microturbulent velocity £ is esti- 
mated from the OH and CO lines. We were able to pro- 
duce satisfactory fits to the data without the inclusion 
of macroturbulence. This is to be expected for obser- 
vations of RSGs with initial masses ~10-20Mq, whose 
macroturbulent velocities are often found to be compa- 
rabl e to the spectral res olution of our observations (see 
also IDavies et all I2008H ) . 

Abundance estimates are obtained by adjusting rel- 
ative abundances in the model atmospheres, comput- 
ing the synthetic spectra, and minimizing the residu- 
als between the observed and synthetic spectra in the 
equivalent widths of certain diagnostic lines. We con- 
centrate on fitting those lines which are relatively un- 
blended and for which reliable atomic/molecular data 
exist. The ro-vibrational lines of OH and CO are used to 
determine oxygen and carbon abundances, while abun- 
dances of other metals are derived from the atomic lines 
of Fe I, Si I, Mg I, Ca I, and Ti I. While the molecule 
CN is included in the analysis, severe blending of the 
lines at the observed spectral resolution prevents an ac- 
curate determination of the N abundance. We do how- 
ever find results that are consistent with CNO equilib- 
rium, that is with N being enhanced by the same level 
that C is depleted. For compilation s of atomic/molecular 
data w e use the Kurucz dat a base 1 , iBiemont fe Grevesse 
1973[) . iHuber fe Herzbergj (fl979l). ISneden fe Lambert 



19821) . and lMelendez fe Barbuvl (l999h . 



For each best-fitting model, we estimate the uncertain- 
ties and solution-robustness by generating a further four 
models with abundances ±0.1dex and ±0.2dex. In addi- 

1 http: //cfa-www. harvard. edu/amdata/ampdata/kurucz23/sekur 



tion, test models are generated in which the stellar pa- 
rameters are varied and the abundances re-tuned, to as- 
sess the degeneracy of the model solutions. In studying 
the residuals between the test-models and the observed 
spectra, we find that solutions of similar statistical signif- 
icance are found with AT c ff±200K, Alogg = ±0.3, and 
A£ = =F0.5km s — 1 , which we take to be the uncertain- 
ties on the stellar parameters. We find that the measured 
abundances are stable across the test models to within 
0.1-0.2dex, while we found that gravities log g = 0.0 ±0.3 
and microturbulcncc £ = 3 ± 0.5 km s _1 produced satis- 
factory fits for all stars modelled. The quoted errors in- 
clude the uncertainty on the continuum placement when 
measuring the line equivalent widths, which can be prob- 
lematic in H-band spectra of cool stars due to the num- 
ber of spectral features. As a further cross-check on our 
results, we re-ran the spectral sysnthesis using the CN 
line-list provided by B. Gustafsson (priv. comm.), find- 
ing negligible differences in derived abundances (see Fig. 

©■ 

Finally, as further validation of our method we note 
that in our recent study of two RSGs in the Galactic Cen- 
tre, we found that our results were consistent with those 
of complimentary studie s of the same object s to within 
our quoted uncertainties (|Davies et al.ll2008bh . while m a 
study of Arcturus we fo und abundances cons istent with 
previous measurements ([Rich fe Or iglia 200§). 

3. RESULTS 

Selected regions of the H-band spectra of one star in 
each cluster are shown in Fig. [TJ Each show a dense ab- 
sorption line spectrum, dominated by neutral lines of Fe, 
several a elements (Si, Ca, Mg, Ti), as well as molecular 
bands of OH and CO. Below, we describe the abundance 
patterns of the element s studied. As r eferen ce, we use 
the solar abundances of lAsplund et alJ (|2005f ). 

The elemental abundances we derive for each star are 
listed in full in Appendix A. To within the uncertain- 
ties, the elemental abundances for each star studied in 
each cluster are consistent. In Table [IJ we show the av- 
erage abundance for each element in each cluster. Here, 
it can be seen that the average abundances in the two 
clusters are similar to within the errors. This evidence, 
together with the proximity of the clusters to one an- 
other (d sep = O.S^Q ^kpc, D08) and their similar ages 
(D07, D08), suggests that the two clusters formed out of 
the same giant molecular cloud. 

We find that, in general, both clusters are slightly 
metal-poor with respect to the solar value. The Fe con- 
tent averaged over the two clusters is [Fe/H] = —0.15 ± 
0.07. The a-elements (Ca, O, Mg, Si, Ti) also show 
abundances which are ~0.1-0.2dex below the solar level. 
These values do not take into a ccount self-de p letion of H 
at the surfaces of the stars: in IDavies et alJ (|2008bf ) we 
show, with the aid of evolutionary models, that the at- 
mospheres of RSGs can be depleted in H by up to O.ldex, 
as the deep convective layers draw the nuclear-processed 
material to the surface. Therefore, we conclude that the 
initial [Fe/H] and [a/H] ratios of the two clusters were 
«0.2-0.3dex below solar. 

The ratio of the a-elements to Fe, which is a powerful 
probe of star-formation history (see Sect. [I]), is consis- 
tent with solar for each of Si, Mg, Ca and Ti . We find 
htmpO/Fc] = ±0.09±0.11 using the lAsplund et all (|2005l ) so- 




wavelength (jum) 

Fig. 1. — : Examples of the H-band spectra of one star in each cluster, showing select regions and the identified lines. 
Observations are shown as dotted lines, our best-fit synthetic spectra are shown as solid black lines. In addition, we also plot 
the computed best-fit spectrum using the Uppsala CN line- list (B. Gustavsson, priv. comm. - solid red lines), with negligible 
differences in the strengths of the diagnostic lines. 



TABLE 1: 

Average elemental abundances of the two clusters, w here A(X) = loe(X/H ) + 12. Columns show (1) the elements 
studied; (2) the mean solar abundances, taken from I Asplund et al.l ( 20051 ) ; (3) & (4) the average abundances for 



RSGC1 and RSGC2 respectively; and (5) the ratio of each element's average abundance from the two clusters 
compared to Fe, normalized to the solar value. The uncertainties in the average cluster abundances are taken as 



the rms standard deviation in each value. 


Element X 


A(X) Q 


^POrsgci 


^P0rSGC2 


[X/Fe] RS GCi 


[X/Fe] RS GC2 


Fe 


7.45 


7.33 


7.28 








± 0.05 


± 0.05 


± 0.03 









8.66 


8.61 


8.59 


0.08 


0.10 




± 0.05 


± 0.07 


± 0.05 


± 0.11 


± 0.09 


Si 


7.51 


7.30 


7.36 


-0.08 


0.02 




± 0.04 


± 0.08 


± 0.08 


± 0.12 


± 0.11 


Mg 


7.53 


7.41 


7.36 


0.00 


-0.00 




± 0.09 


± 0.13 


± 0.09 


± 0.17 


± 0.14 


Ca 


6.31 


6.20 


6.18 


0.02 


0.04 




± 0.04 


± 0.05 


± 0.04 


± 0.10 


± 0.08 


Ti 


4.90 


4.92 


4.92 


0.14 


0.19 




± 0.06 


± 0.10 


± 0.06 


± 0.14 


± 0.08 


C 


8.39 


7.95 


7.89 


-0.32 


-0.33 




± 0.05 


± 0.07 


± 0.06 


± 0.11 


± 0.10 



The metallicity of the Scutum RSG Clusters 



5 



lar abundances. When the iGrevesse fe Sauval (|1998l ) so- 
lar values are used, we find [O/Fe] = — 0.03±0.11. Given 
the recent fluctuations in the derived solar O abundance, 
at this time we conclude for now that the O/Fe ratio in 
the clusters is consistent with solar. 

The levels of carbon in the stars' atmospheres appear 
to be significantly depleted, consistent with their evolved 
status. This result is expanded upon in detail in Sect.[4j 

4. COMPARISON WITH STELLAR EVOLUTIONARY 
MODELS 

The results presented in Sect. [3] allow us to test vari- 
ous predictions of stellar evolutionary models. As stars 
enter the RSG phase, the products of CNO burning are 
expected to be seen at the surface, specifically the en- 
richment of N at the expense of C. From the values in 
Table [3j it can be seen that all stars in both clusters 
show evidence of C depiction. In Fig. [2] we investigate 
quantitatively the predictions for surface C abundances 
as a function of evolutionary state. 

Throughout our analysis, we consider the C fraction 
in relation to the Fe abundance. In evolutionary models, 
the relative abundance of an element X with respect to 
hydrogen X/H depends on two input parameters: the 
relative abundances of the heavy elements, and the ratio 
of the heavy element abundances to H (i.e. the metal- 
licity, Z). However, the ratio X /Fe depends only on 
the relative abundances of the heavy elements. Thus, by 
comparing C as a function of Fe, we reduce the model- 
dependency of our analysis. As the relative abundances 
of the heavier elements in the stars studied here appear 
to be consistent with solar (see Table QJ , we use models 
with solar-scaled heavy element abundances. 

4.1. Non-rotating models vs rotating models 

In Fig. [5] (left) we plot the C/Fe fraction as a function 
of luminosity for the stars observed. The stellar lumi- 
nosities are taken from D07 and D08. Overplotted are 
mass tracks at initial masses of 12M Q , 15M Q and 2OAf 
from the non-rotating Geneva models. For reference, the 
initial masses of the RSGs in each cluster (derived from 
the rotating Geneva models) are 18M Q for RSGC1 and 
UM Q for RSGC2 (D08). As we find a sub-solar iron 
content for the stars in each cluste r, we plot tracks us- 
ing both solar metallicity (Z=0.02 ISchaller et al1l!992f) 
and Z=0.008 ()Schaerer et al.l I1993D . The iron content 
we derive for the RSG clusters indicates a metallicity 
somewhere between these two values. In the right-hand 
panel of Fig. [21 we plot the same quantities in com- 
parison to the Geneva models which include rotation , 
at Z = 0.02 and Z = 0.004 (|Mevnet fe Maederl l2000l : 
iMaeder fe Mevnedl20"ol respectively) 2 

When comparing our derived C abundances with those 
predicted by the Geneva evolutionary models, it is impor- 
tant to first note that the se models use the relative h eavy 
elemental abundances o f IGrevesse fe Sauvall (fl998h . In 
the 3-D solar model of lAsplund et al.1 1 20051 ) the frac- 
tional abundances of C and Fe have been revised slightly. 
For this analysis we assume that this small change in rel- 
ative metal abundances would have no major impact on 
the output of the Geneva evolutionary code other than 

2 Geneva rotating models at Z=0.008 are unavailable at the time 
of writing. 



to cause a linear displacement in abundance-space. We 
have shifted each mass track by 0.105dex- the magnitude 
of the revision in the Solar [C/Fe] fraction. We indicate 
the magnitude of this displacement with an arrow in the 
upper-right of each of the panels in Fig. [21 

In general, the non-rotating models do not quite re- 
produce the observed level of C depletion in the RSGs 
of the two Scutum clusters, even when the updated solar 
[C/Fe] is taken into account (Fig. [21 left). There seems to 
be much better agreement between the data and the ro- 
tating models, which are more heavily depleted than the 
non-rotating models due to rotationally-enhanced inter- 
nal mixing. Sub-Solar metallicity models suffer a slightly 
greater level of C depletion, due to decreased mass-loss 
on the main sequence which inhibits the loss of angu- 
lar momentum and increases the level of internal mixing. 
This can be seen in the right-hand panel of Fig. [21 where 
the low-Z mass-tracks have systematically lower C /Fe ra- 
tios. While we are unable to discern between the Solar 
and sub-Solar mass-tracks, we can say that the rotating 
models provide a much better fit to the [C/Fe] levels, 
while it has already been shown that these models do a 
better job of reproducing the luminosity distributions of 
the two clusters (D08). 

4.2. Binary models 

To assess the impact of binarity on surface abundances 
we refer t o the Cambridg e binar y star evolution code pre- 
sented in E ldridge et al.1 (|2008l ). The code includes the 
effect of Roche Lobe Overflow and Common-Envelope 
Evolution, processes which are particularly important 
when the primary star evolves to the RSG phase due 
to the expanded atmosphere. 

Briefly, the evolution to the RSG phase of a massive 
star can result in the star filling its Roche Lobe and trans- 
ferring matter to the companion star. This reveals the 
lower layers of the primary's atmosphere in which the ef- 
fects of CNO burning are more evident. However, mass- 
transfer of the primary's envelope onto the secondary 
also restricts the primary's redward evol ution, resulting 
in dra stically shortened RSG lifetimes. Eldridg e et al.l 
(2008) showed that, for a sample of binary systems with 
broad distributions of separations and mass-ratios, the 
mean RSG lifetime averaged across the sample was re- 
duced by a factor of ~3. However, for low-separation bi- 
nary systems the RSG lifetime-reduction is clearly much 
more severe. Therefore, the chances of observing a RSG 
in a close binary are extremely small, given the brevity 
of the phase in single stars (ijtsG ~lMyr for a single star 
in the Eldridge et al. models). 

While binary models are able to increase the level of 
C depletion, the drastically reduced RSG lifetimes lead 
us to conclude that close binary interaction cannot have 
played a significant role in the evolution of the RSGs in 
the two clusters. This is not to say that binarity is not 
prevalent in the clusters, more that the RSGs in the two 
clusters that we see today are unlikely to have evolved in 
a close binary system. This has two important implica- 
tions, which we now discuss. 

Firstly, in F06, D07 and D08, the masses of the clus- 
ters were estimated by counting the numbers of RSGs, 
comparing to population synthesis models created using 
single star evolutionary codes, and extrapolating over the 
rest of the Initial Mass Function. If the clusters have 
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Fig. 2. — : Carbon fraction as a function of stellar luminosity. Overplotted are mass-tracks from the Geneva non-rotating (left) 
and rotating (right) isochrones at solar and sub-solar metallicities (see text for references). Note that the sub-solar metallicities 
are different in each panel. The mass-tracks have been shifted downwards to reflect the recent revision in the Solar [C/Fe] - the 
arrow in the top-right of each panel indicates the magnitude of this change (see text for details). 



high binary fractions, then a significant number of cluster 
members which had the same initial masses of the RSGs 
we see now will have skipped the RSG phase. The not- 
unreasonable parameters for a binary distribution used 
by Eldridge et al. suggest a reduction in the number of 
RSGs by a factor of ^3 compared with single-star mod- 
els. Therefore, the total cluster masses may have been 
underestimed by a similar factor, implying initial masses 
in excess of 10 5 M Q , which would make them the most 
massive young clusters in the Galaxy. 

Secondly, the mere fact that we are looking at RSGs 
in these clusters suggests that we are dealing with stars 
which have in effect evolved in isolation. That is, it is rea- 
sonably safe to assume that the stars are single. This is 
an extremely important consideration if we are to use the 
stars to measure the velocity dispersion of the clusters, 
and hence determine their virial masses. Contamination 
by the high orbital velocities in close massive binaries 
would lead to a significant systematic overestimate of a 
sample's velocity dispersion, and hence of the virial mass. 
On the one hand, this should serve as a note of caution 
when interpreting the velocity dispersions in young mas- 
sive clusters when measured from blue stars, particularly 
as the massive binary fraction of W esterlund 1 has re - 
cently been inferred to be very high (|Clark et al.ll2008f ). 
On the other hand, this would also seem to suggest that 
RSG spectral features provide the most reli able measure- 
ment of cluster velocity dispersions (see e.g. lMengel et al.l 
120081) . 

5. THE RSG CLUSTERS AND THE STAR-FORMATION 
HISTORY OF THE GALAXY 

In this section, we compare the relative abundances de- 
rived for the RSG clusters with the standard indicators 
of star- formation history, namely the ratios [Fe/H] and 
[a/Fe]. In broad terms, the primary source of a-element 
enrichment in the Galaxy is through Type-II supernovae 
(SNe), that is the core-collapse of stars with masses 
greater than ~8M Q . On the other hand, the main source 
of Fe enrichment is SNe la - thermonuclear explosions 
of low mass stars which reach the Chandrasekhar limit 
through accretion from a companion. Hence, as Fe- and 



a-enrichment occurs on vastly different timescales, star- 
formation histories can be probed through comparisons 
of the a and Fe abundances 

5.1. The abundance patterns of the Galaxy 

In the Galactic 'thin' disk, with Galacto-centric dis- 
tances Rqc ^4kpc and scale height \h\ ^ lOOpc, it 
has been found from abundance analysis of Cepheids 
and late- type dwarfs that the ratio [a/Fe] is app roxi- 
mately solar (|Bensbv et alJl2004t iLuck et al.ll2006l and 
references therein). Meanwhile, in the 'thick' disk 
(\h\ p> lOOpc), halo and bulge, [a/Fe] is found to be 
super-solar from analyses of Red Giants (iRich fe Origlial 
2001 iCunha fc Smi th 2006; Fulbrigh tet all 120061 120071 



Lecureur et alJl2007t iMelendez et alJl2008f > 



In addition, many authors have found evidence for 
a trend of increasing chemical abundances with de- 
creasing Rgc, the so-called Ga lactic metallicity gra- 



dient. Ana l yses of B dwarfs (iRolleston et al. 120001 : 
iSmartt et al]l200lT). Cepheids (ILuck et al.l 120061). UC- 



HII regions ( e.g. [A fficrbach et al. 19971) . and Planetary 
Nebulae (e.g. iMaciel fe Quirezal ll999) have all found ev- 
idence for increasing metal abundances, both of Fe and 
a-elements, at lower Rgc- Further, the res ults of An- 
drieysky and collaborators (summarised in ILuck et al.l 
2006), who studied the abundances of both Fe and the 
a-element Ca in Cepheids, showed that [Ca/Fe] was so- 
lar and roughly constant throughout the thin disk out to 
Rgc ~ 8kpc, with tentative evidence for an increase in 
[Ca/Fe] in the outer Galaxy. 

These abundance patterns are thought to be a result 
of the Galaxy's evolution and star-formation history: 
the bulge and hal o formed very quickly (~0.1-0.5Gyr, 
iBallero et al.l l2007) , during which time the chemical evo- 
lution was dominated by SNe II. After this time, star for- 
mation ceased and the a-heavy abundances were frozen- 
in. However, in the Galactic disk star-formation contin- 
ued at some rate which was a function of Rgc- Even- 
tually the enrichment of the Fe abundance by SNe la 
brought the a/Fe ratio to below that of the halo and 
bulge. 
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Fig. 3. — : The trend of Fe/H abundance in the inner 6kpc 
of the Galaxy. The dot ted line in each panel indicates the rel- 
evant solar value from lAsplund et al.l J2005). The references 
for the literature data are given in Tabled 
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Fig. 4. — : Same as Fig. [3] but for oxygen (top panel), 
magnesium (middle panel) and silicon (bottom panel). 

5.2. Comparison between Galactic and RSGC 
abundances 

5.2.1. Iron 

As stated in Sect. 02 we find that the average Fe/H 
abundance ratios for the two clusters are 0.1-0.2dex be- 
low solar. This is neglecting the effects of self-depletion 
in the atmopsheres in RSGs, which we estimate to be a 
further ~ —0.05 — O.lde x for stars with initi al masses of 
14-18M© (see Fig. 2 of IDavies et al]|2008bf ). 



TABLE 2: 

References for the data shown in Figures [3] - [6] 



Object 



References 



B stars 
Hll regions 
Other RSGs 
LBVs 
Cepheids 



Smartt et al. (2001); Dafion & Cunha (2004) 

Afflcrbac h et aT7T997T ~ 

Davies et al. f 2008bT~Cunha et al. (2007) 

Naiarro et al (2008) " 

Andrievskv et al. (20021 



In Fig. [3] we compare the Fe/H ratios of the 
clusters with other Fe/H measurements in the inner 
Galaxy. In anal ysis of the atmosp heres of Cepheids, 
lAndrievskv et al.l (|2002f ): lLuck et all (|2006t ) find that the 
Fe/H ratio increases towards the Galactic Centre (GC) at 
a rate of -0.06dexkpc _1 , reaching «+0.2dex above solar 
at Rqc ^ 6kpc. In the GC, recent studies suggest that 
the Fe abundance i s roughly solar (iNaiarro et al.l l2008t 
IDavies et all l2008bl e.g.). We note that ICunha et all 
(12007J) found slightly super-solar Fe levels from studies 
of RSGs in the GC. However, the initial Fe/H ratios of 
Cunha et al.'s sample were likely 0.05-0.1dex lower than 
reported due to the affore-mentioned self-depletion effect. 

From Fig. [3] it is clear that the RSGCs do not fol- 
low the radial Galactic Fe/H trend - the clusters ex- 
hibit Fe abundances which are ~0.5dex below those of 
the Cepheids in the inner disk. This plot may simply 
illustrate the inadequacies of a one-dimensional param- 
eterization of abudance trends, particularly in the inner 
Galaxy. The presence of physical structures in the in- 
ner regions of the disk, such as the Bar, and the ends of 
the spiral arms, represent large azimuthal fluctuations in 
stellar/gas density and star- formation rates. As such, it 
seems unlikely that chemical abundances should be ho- 
mogeneous around the inner disk. This is a point that 
we investigate further in Sect. 15.231 



5.2.2. The a- elements: O, Ca, Si, Mg, and Ti 

As with Fe, these elements are again found to be sub- 
solar in relation to H. However, the ratio of a/Fe is 
solar to within the errors, consistent with a disk pop- 
ulation. In Fig. |4] we again put these results in con- 
text with measurements of objects at similar Rgc- In a 
search of the literature, we have found abundance mea- 
surements in the inner di sk for O, Mg and Si for Cepheids 
(lAndrievskv et al.ll2002ft . B dwarfs (jSmartt et all 120011 : 
iDaflon fe Cmihan 2004) ; while we also find O measure- 
ments from H II regions (|Afflerbach et al.l [19971) . The 
H II region measurements are less reliable as they are 
highly model-dependent - they determine temperatures 
by matching the amount of ionizing UV flux (inferred 
from the radio) to model stellar atmospheres, which in 
itself depends on metallicity. Nevertheless, we include 
their data for completeness. 

Figure [4] again shows that the RSG clusters appear 
to have systematically lower abundances than other ob- 
jects located at similar Galacto-centric distances. How- 
ever, when the B dwarf measurements are included, the 
discrepancy of the RSGC data appears less striking. In- 
deed, taken together, the cepheid, RSGC and B dwarf 
measurements seem to suggest that there be large abun- 
dance variations in the inner Galaxy. 
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Fig. 5. — : Location in the Galaxy of the B stars, H n regions 
and Cepheids used as probes of the metallicity gradient. We 
assume the Galactic centre to be 8kpc from Earth. The dot- 
ted circles show Galacto-centric distances of 2, 4, and 6kpc. 
The references for the literature data are given in Table [2] 

5.2.3. Azimuthal variations in relative abundances 

As stated in the preceeding sections, the simple one- 
dimensional parameterization of abundance patterns in a 
metallicity gradient may be too obtuse when dealing with 
the inner regions of the Galaxy. Indeed, the presence of 
the Bar at Rgc & 4kpc implies large azimuthal fluctua- 
tions in stellar/gas density, hence it is not inconceivable 
that variations in abundances may also exist. 

In Fig. [5J we show the locations of the objects plot- 
ted in Figs. [3] and 2] in the X-Y plane. We see that the 
Cepheids and the RSGCs, which have abundance dif- 
ferences of ~0.4dex, are separated in physical space by 
approximately 4kpc. Curiously, the B dwarfs and H II 
regions, which have intermediate chemical abundances 
between the RSGCs and Cepheids, are also located be- 
tween the two in physical space. This raises the intrigu- 
ing possibility of large-scale (~kpc) azimuthal abundance 
variations in the inner disk. 

To explore this prospect further, in Fig. [5] we plot the 
relative abundances of the Cepheids, B dwarfs, Hn re- 
gions and RSGCs as a function of their azimuthal angle 
about the Galactic Centre, <^gc- We choose to plot the 
abundances as a function of the cosine of the angle, as 
any periodic variations should show up as linear trends in 
cos(r/>Gc)-space. We define the positive direction as being 
clockwise when viewed from above, and we have rotated 
the coordinate system such that the RSG clusters, which 
have the lowest abundances, have cos ((/>') = 1, where 
4>' = 4>gc ~ 0R.SGC- In practice this coordinate system 
transformation made very little difference to the results. 

From Fig. [5] we see that, though there is no strong ev- 
idence of a gradient in Mg, trends do seem to exist in O 
and Si. This is true even if the model-dependent H n- 
region points are removed. We find that the gradients in 
O and Si, determined from linear regressions, are signifi- 
cant at the ~5cr level. Further, the gradients are similar 
in each element to within the uncertainties. When the 
coordinate system was not rotated (see above) , this gra- 
dient significance was still in excess of 4cr. 

These measurements are probing a limited range in 
azimuth, so these observed gradients are representative 
of localised fluctuations in abundance levels in the in- 
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Fig. 6. — : Azimuthal abundance gradients for oxygen, 
magnesium and silicon, created using objects with Galacto- 
centric distances of 3.5 to 6.0 kpc. The coordinate system has 
been rotated such that (4>gc — 0rsgc) = 0°. The references 
for the literature data are given in Table [2] 

ner disk. Moreover, these fluctuations apparently coin- 
cide with the end of the Bar. Given that we find these 
statistically-significant variations, we identify three sta- 
tistical explanations for the observed trends which we 
discuss below. 

Artifacts of random sampling If the abundances in the 

inner Galaxy are uniform, but have large scatter, then a 
gradient in abundances may be artificially inferred if we 
happen to sample the extremes of the abundance distri- 
bution at the extremes of the observed cpac range. This 
would seem like a plausible explanation for the measured 
gradient in 0gc, as the data-points in Fig. [5] are few and 
poorly sampled at the extremes of 4>gc- 

To investigate quantitatively the likelihood of this sce- 
nario, we have performed monte-carlo simulations to gen- 
erate random abundance levels at the same sampling in 
cos(</>')-space as the data plotted in Fig. [6] For each 
simulation, we generated data with uniformly random 
distribution with a peak-to-peak amplitude of 0.85dex 
- similar to that of the observed data. We set the un- 
certainty in each value of A(X) to the same as in the 
observations at that value of cos(</>') 3 . We calculated the 
gradient g and uncertainty o~g of the simulated data, and 
measured the significance of the gradient Q = g/ug. The 

3 Setting all measurement uncertainties to ±0.15dex had no ef- 
fect on the results of our analysis. 




Fig. 7. — : The probability of obtaining a gradient with significance Q (= g/crg ) for three different values of AA, which is 
defined as systematic difference in the measured abundance of an element when determined from two different methods. We plot 
simulations for two different intrinsic abundance scatters: AZ — ±0.2dex (left) and AZ = ±0.5dex (right). The significance of 
our measured azimuthal gradients in O and Si, Q — 4, is indicated by the dotted line. Each plot shows that, for larger systematic 
differences, the probability of artificially inferring a gradient is higher. Meanwhile, comparison between the two figures shows 
that gradients are more likely to be inferred if the intrinsic scatter is larger. 



simulation was repeated 10 5 times and the frequency of 
obtaining a measured gradient with significance Q > 4 
was recorded. 

We found that the likelihood of measuring a gradient 
at the Q — 4 significance level in the O and Si data was 
low, ~14%. As it may be argued that the two RSGC 
measurements may not be independent of one another, 
we then repeated the experiment, this time forcing the 
two RSGC data-points to the same value. This time, the 
probabilities of obtaining a Q — 4 gradient rose to «23% 
for each of O and Si. We repeated this analysis, this time 
generating random A(X) values from a guassian distri- 
bution with standard deviation a A = 0.3, the observed 
spread of the Si abundances. This time, we found that 
the probability of obtaining Q = 4 was 22%, or 30% 
when the RSGC points were forced to the same value. 

From this analysis, we conclude that the measured gra- 
dients in Si and O are unlikely to be a result of random 
sampling. That is to say, if large azimuthal abundance 
dispersions exist in the inner disk, the likelihood of arti- 
ficially inferring the gradients we measure from random 
sampling is <;30%. 

Systematic offsets between analysis methods - - Another 
explanation for the observed gradient may be that the 
abundances in the inner disk are uniform, but a system- 
atic offset exists between abundances measured using dif- 
ferent methods. For example, Si abundances measured 
from Cepheids may yield values which are systematically 
higher than those measured from RSGs, even though the 
true abundances of the two objects are the same. Indeed, 
one may argue from the bottom panel of Fig. [H] that the 
Si abundances of the Cepheids are constant to within the 
errors, as are the two RSG clusters, with the two classes 
of objects offset from one another by ^0.5dex. 
We performed a quantitative investigation of this, 



again using monte-carlo simulations. We assumed that 
the abundances Z in the inner disk are uniform, with 
an intrinsic peak-to-peak scatter of AZ. We randomly 
generated abundance values at the same distribution in 
cos(</>') as the data in Fig. [6l We then artificially raised 
the abundances of the cepheid points by AA, and de- 
creased the RSGC points by the same quantity. That is, 
a value of AA = ±0.2 separates the two classes of ob- 
ject by 0.4dex in abundance space. We again forced the 
RSGC abundaces to be the same, as these measurements 
may not be independent of one another. The abundances 
of the B stars were randomized, with a distribution of 
±0.4dex about the mean value, similar to the data in 
Fig. [Sj At the end of each simulation, the gradient g 
and its significance Q was measured, and the experiment 
repeated 10 5 times. 

In Fig. we plot the likelihood of measuring a sta- 
tistically significant gradient for three different values of 
AA, using the same cos(^>') values as those in the bot- 
tom panel of Fig. [6l We then reproduce this analysis 
for two different intrinsic abundance scatters, AZ . For 
example, the right-hand panel of Fig. [Jj shows the frac- 
tion of simulations for which a gradient with significance 
Q is measured for data with a peak-to-peak abundance 
scatter of AZ =0.5dex. In this panel, for a systematic 
offset AA = ±0.1dex - the same as quoted uncertainties 
in the data - the probability of inferring a Q > 4 gra- 
dient is 20%. For an offset of AA = ±0.2, which is the 
value typically quoted as being the absolute uncertainty 
on any abundance measurement, the probability is raised 
to 50-60%. 

From this analysis, we conclude that it is possible that 
the azimuthal abundance gradient we measure is due to 
systematic offsets between differing measurement tech- 
niques. However, for a 50% probability of obtaining 
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Fig. 8. — : Azimuthal [O/H] variations in NGC 4736. Sym- 
bol sizes represent the magnitude of the abundance levels in 
the various H uregions that surround the central region, as in- 
dicated by the legend. The shaded region denotes the rough 
extent and position angle (110±30°) of the galaxy's bar. 

the Q = 4 gradients we observe, the offset between the 
RSGCs and the Cepheids must be large, ~0.4dex (i.e. 
AA = ±0.2). We note that in the Galactic Centre, we 
found good agreement (better than ±0.2dex) between 
our analysis of RSGs and that of ot her classes of objects 
such as Luminous Blue Variables (jNajarro 
iDavies et ai1l2008bD . 

Real trends Having studied two scenarios that could 

lead us to artificially measure a gradient where none ex- 
isted, we now discuss the possibility that the measured 
azimuthal abundance variations are real. If such be- 
haviour exists in our Galaxy, then it is reasonable to 
assume that it should also be present in external galax- 
ies. In the following section, we investigate the presence 
of azimuthal abundance variations in the face-on spiral 
NGC 4736. 

5.2.4. Case-study: NGC 4736 (= M94) 

In looking for an azimuthal abundance gradient in an 
external galaxy, we have chosen this object for several 
reasons. Firstly, while not having the exact same type 
as our own Galaxy, it is a spiral with a central bar, 
and is relatively face-on (|Block et al.lll994l ). Secondly, 
it has a ring of H n regio ns in the inn er disk, which are 
well sampled in azimuth ()Lvnds![l974f ). Thirdly, oxygen 
abundan ces of these Hli r egion s are present in the lit- 
erature: iMartin fc Bellevl (|1997l ) studied the H II region 
ring with narrow-band filters to measure the fluxes of 
certain diagnostic lines in the optical, and used them to 
derive abundances. However, the authors did not use 
the abundances to study azimuthal variations, only the 
presence of a radial gradient. This study utilized a homo- 
geneous methodology, and so we do not have to take into 
account systematic errors induced by differing analysis 
techniques. 

In Fig. [8] (left) we plot the locations of the H n regions 
in the X — Y plane, corrected for the galaxy's inclina- 
tion. The size of the plotting symbols indicates the mag- 
nitude of the abundance of each H n region with respect 
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Fig. 9. — : The abundance ratio [O/H] as a function of az- 
imuthal angle, with the coordinate system rotated such that 
the central bar (0/ rm( , ar =llO±3O°) is oriented at cos(<^>') = 1. 
In the top panel, we plot [O/H] versus azimuthal angle in the 
clockwise direction away from the bar, plotting those points 
with (f> >80 c '(=0bar — c<^bar)- Similarly, in the bottom panel 
we [O/H] as a function of cos(0) in the anti-clockwise direc- 
tion, plotting those points with <f> < 140° (= (<^b ar + 0"</>bar) 

to the global average. To the eye, it seems that the abun- 
dances in the lower-right quadrant of the figure tend to 
be lower than average; while at the left edge of the ring 
at (X,Y) w (—1,0) the abundances are highest. Inter- 
estingly, the minimum in abundance levels corresponds 
roughly with the end of the bar (</>bar = 110±30° indi- 
cated by the shaded region in Fig. [8]- iBlock et alJll994t 
IWong fc Blit3l2000D . This is similar to the trend seen in 
our own Galaxy, where the abundance minimum is found 
in the RSG clusters located at the end of the Galaxy's 
Bar. 

In Fig. [9l we plot the abundance levels as a function 
of the cosine of the azimuthal angle. Again, as in Fig. [6] 
we have rotated the coordinate system to such that the 
minimum abundances correspond to cos(0) = 1.0. For 
this rotation angle we adopt the position angle of the bar, 
110±30°. The top panel of Fig. [9] shows the abundance 
trend as a function of cos(^) in the clockwise direction 
from the lower limit of the bar's orientation (110°- 30°= 
80°). Similarly, in the bottom panel we show the abun- 
dance trend in the anti-clockwise direction up to the up- 
per limit of the bar's orientation (140°). We opt to plot 
the abundance variations in the two directions separately, 
as the trends in each direction are not necessarily similar. 

In each panel of Fig. [5] we make a linear fit to the [O /H] 
distribution and indicate on the plot the measured gra- 
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dient and its uncertainty. The maximum position angle 
that we have chosen to plot in each panel, <^6 or ±160°, was 
somewhat arbitrarily chosen to illustrate a wide range of 
azimuth as well as include a sufficiently large number of 
data-points, though we note that the range chosen does 
not significantly affect Q-values of the measured gradi- 
ents. 

From the two panels of Fig. [9] we see that abun- 
dance gradients are measured at the Q — 4 (clock- 
wise) and Q = 2 (anti-clockwise) levels. Though ten- 
tative, this does provide some support for the idea of lo- 
calised azimuthal abundance variations in the inner disk 
of NGC 4736. We again state that, as in our analysis of 
the inner Galaxy, these 'gradients' are more likely to be 
indicative of systematic undulations in abundance lev- 
els rather than true sinusoidal variations with azimuthal 
angle. 

5.2.5. Bars, and azimuthal abundance variations 

To summarize the results of the previous section, we 
find it unlikely that the azimuthal abundance patterns 
we observe in the inner Galaxy are an artifact of random 
sampling. We cannot rule out that the observed vari- 
ations result from systematic offsets between different 
methodologies. However, in using data from a homoge- 
neous study of the external galaxy NGC 4736 where we 
are able to make a differential analysis, we find tentative 
evidence for similar abundance variations. In each case, 
the lowest abundances were observed at the end of the 
central bar. 

Given these apparent azimuthal variations in abun- 
dances seen in both our Galaxy and NGC 4736, we 
now discuss how such could patterns arise. It is known 
that there is a link between abundance trends and the 
presence of bars - radial abundance gradients tend to 
be flatter in barred spirals than in un barred spirals 
(|Zaritskv et alj|1994t iMartin fc Rovlll994f) . Indeed, the 
observed abundance gradient in our Galaxy is typical of 
barred spirals, and much shallower than unbarred spi- 
rals. The common explanation for this is that radial gas 
motions induced by the bar potential smooth-out any 
pre-existing abundance gradient on timescales of a few 
x lOOMyr. However, this can introduce significant az- 
imuthal variations, especially between the spiral arms 
and the int er-arm regions, due to intense but patchy star- 
formation (jFriedli et al.lll994 iFriedli fc Benzlll995l ). 

In NGC 4736, kinematic studies have shown that the 
H II ring is coincident with the outer Lindblad resonance 
of the bar and the inner Lindblad resonance of the central 
bulge. The ring itself is likely a result of radial inward 
gas motions ou tside of this resonance zone and outward 
motions inside (jWong fc Blitzll2000l ; iMunoz-Tunon et all 
l200l . The ensuing star-formation within the ring, as 
the colliding material is shocked, has been shown to 
be pa tchy from CO and Ha observations (jWong fe Blitz! 
2000) . As the timescale for chemical enrichment by mas- 
sive stars is comparable to the dynamical timescales in 
the region of the bar (few xlOOMyr), it is conceivable 
that notable fluctuations in the azimuthal abundances of 
the a-elements would survive. Similarly, in our Galaxy, 
azimuthal abundance variations could be produced by 
intense but patchy star-formation in the region of the 
Bar's outer Lindblad resonance at i?QC~4kpc. Indeed, 
this starburst episode could be responsible for creating 



the RSG clusters. However, it is unclear how this effect 
would show itself in the Fe abundances, as the enrich- 
ment timescale for this element is longer (~ 10 8 — 10 9 yr). 
Clearly, the abundances in the inner Galactic disk re- 
quire further study, and would benefit from a homoge- 
nous s tudy such as that of NGC 4736 bv lMartin fc Belleyl 
(|1997f ). Massive young clusters such as RSGC1 and 
RSGC2, from which abundances of several elements may 
be determined in the near-IR, may represent the best tool 
with which to do this. Though known examples of such 
clusters are currently low in number, searches of infrared 
Galactic plane surveys have revealed many more candi- 
dates, and the number of known young massive clusters 
in the inner Galaxy looks set to rise in the next few years. 

6. CONCLUSIONS 

We have used high-resolution near-infrared spec- 
troscopy and detailed spectral synthesis to determine the 
chemical abundances of the two Scutum Red Supergiant 
clusters, RSGC1 and RSGC2, which are located in the 
inner Galaxy close to the end of the Bar. Our results can 
be summarized as follows: 

• The stars atmopsheres are carbon-depleted, and 
the average levels of depletion agree well with the 
predictions of rotating single-star models. Non- 
rotating models predict C/Fe surface abundances 
which are higher than observed. 

• The average iron abundances (Fe/H) are found to 
be sub-solar by 0.2-0.3dex. This is in apparent con- 
flict with the Fe abundances of objects at similar 
Galactocentric distances, which have typically been 
found to be ^0.5dex higher. 

• The abundances of the a-elements are also sub- 
solar. The ratio of a/Fe is fully consistent with 
solar, as is typical for objects in the thin disk. 

To explain the apparently unusual abundances of the 
RSG clusters with respect to other objects in the inner 
disk, we have collected several abundance measurements 
of different classes of object in the literature. It appears 
that large azimuthal abundance variations exist in the 
inner disk, as is predicted by theory for the secular evo- 
lution of galaxies with bars. Further, we find that the 
objects in the inner disk with strongly super-solar abun- 
dances are typically located at negative Galactic longi- 
tude /, while the RSG clusters are located at positive 
I. We find tentative evidence for large-scale azimuthal 
abundance variations in the data, significant at the 4- 
5<7-level. The presence of these variations can only be 
explained away if abundance studies of different types of 
object have significant systematic offsets from one an- 
other. In a differential analysis study of the external 
barred galaxy NGC 4736 using homogeneous methodol- 
ogy, we again find tentative evidence for large-scale az- 
imuthal abundance fluctuations. We suggest that these 
results for both our Galaxy and NGC 4736 are consistent 
with the predictions of the Friedi & Benz model for chem- 
ical evolution in barred spirals, whereby the bar potential 
which causes a pile-up of material at the outer Lindblad 
resonance results in intense but patchy star-formation in 
the inner region of the disk. As the dynamical timescale 
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is comparable to the timescale for chemical enrichment, 
large azimuthal abundance variations can form. 
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Ap pendix: The che mica l abundances of the stars observed in each cluster. Star identification numbers are taken 
from lFiger et alj (|2006f) and lDavies et all (|2007f ) 
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TABLE 3: 

Elemental abundances for each star studied. Abundances for each element 
X are expressed in the form A(X) — log(A/H) + 12. For reference, the 
solar abundances derived by Asplund et al. (2005) are also shown. 



Star 


T c ff 


A(Fc) 


A(O) 


A(Si) 








A (C\ 


Solar Abundances (Asplund et al. 2005) 














7.45 


8.66 


7.51 


/ .Oo 


a q 1 


a on 


o.oy 






± 0.05 


± 0.05 


± 0.04 


4- n no 
it u.uy 


4- O O/l 
It U.U4 


4- n nfi 
± u.uo 


4- n r\K 
it U.UO 


RSGC1 Abundances 














1 


3600 


7.36 


8.58 


7.21 


7.57 


6.16 


4.98 


7.82 


— 




± 0.09 


± 0.08 


± 0.14 


± 0.12 


± 0.07 


± 0.10 


± 0.09 


2 




7.30 


8.59 


7.25 


7.28 


6.19 


4.72 


7.92 






± 0.08 


± 0.09 


± 0.14 


± 0.12 


± 0.07 


± 0.10 


± 0.09 


3 


■ 


7.40 


8.58 


7.30 


7.48 


6.24 


5.02 


7.92 






± 0.09 


± 0.13 


± 0.15 


± 0.13 


± 0.08 


± 0.12 


± 0.09 


A 
1 


qonn 


7 IK 
1 .OO 


8.76 


7.32 


7.36 


6.28 


5.02 


7.92 






± 0.09 


± 0.07 


± 0.11 


± 0.12 


± 0.07 


± 0.08 


± 0.09 


5 


3600 


7.29 


8.63 


7.22 


7.28 


6.20 


4.92 


7.97 






± 0.08 


± 0.10 


± 0.14 


± 0.12 


± 0.07 


± 0.10 


± 0.09 


6 


3400 


7.41 


8.58 


7.35 


7.64 


6.29 


5.02 


8.02 






± 0.08 


± 0.12 


± 0.15 


± 0.13 


± 0.08 


± 0.12 


± 0.09 


7 


3600 


7.27 


8.54 


7.25 


7.28 


6.12 


4.82 


7.92 






± 0.08 


± 0.09 


± 0.14 


± 0.12 


± 0.07 


± 0.09 


± 0.09 


8 


3600 


7.32 


8.61 


7.25 


7.27 


6.18 


4.92 


8.02 






± 0.08 


± 0.09 


± 0.14 


± 0.12 


± 0.07 


± 0.10 


± 0.09 


9 


3400 


7.37 


8.53 


7.25 


7.54 


6.19 


4.92 


7.92 






± 0.09 


± 0.11 


± 0.15 


± 0.13 


± 0.08 


± 0.10 


± 0.09 


10 


3600 


7.31 


8.72 


7.35 


7.48 


6.22 


5.02 


8.02 






± 0.09 


± 0.10 


± 0.14 


± 0.12 


± 0.08 


± 0.11 


± 0.09 


11 


3600 


7.29 


8.57 


7.35 


7.38 


6.21 


4.92 


8.02 






± 0.08 


± 0.10 


± 0.14 


± 0.12 


± 0.07 


± 0.10 


± 0.09 


12 


3800 


7.38 


8.70 


7.53 


7.53 


6.26 


4.92 


8.02 






± 0.09 


± 0.07 


± 0.11 


± 0.12 


± 0.06 


± 0.07 


± 0.09 


13 


3800 


7.27 


8.68 


7.25 


7.28 


6.18 


4.72 


7.82 






± 0.07 


± 0.07 


± 0.11 
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± 0.09 


14 


3600 


7.25 


8.54 


7.35 


7.38 


6.13 


4.92 


7.92 






± 0.08 


± 0.10 


± 0.14 


± 0.12 


± 0.07 


± 0.10 


± 0.09 



RSGC2 Abundances 



2 


3600 


7.04 


8.37 


7.05 


7.08 


6.01 


4.73 


7.72 






± 0.08 


± 0.09 


± 0.13 
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± 0.06 
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± 0.09 
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3600 


7.22 
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6.11 


4.92 


7.92 






± 0.08 
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± 0.07 
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± 0.09 


5 


3600 


7.29 


8.54 


7.25 


7.25 


6.18 


4.92 


7.92 






± 0.08 


± 0.10 


± 0.14 


± 0.12 


± 0.07 


± 0.10 


± 0.09 


6 


3600 


7.30 


8.60 


7.50 


7.53 


6.21 


4.92 


7.82 






± 0.08 
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± 0.14 


± 0.12 


± 0.07 
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± 0.09 


8 


3600 
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6.23 


4.92 


8.02 
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3600 
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10 
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14 
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7.26 
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